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ABSTRACT 

The Cosmic Evolution Survey (COSMOS) allows for the first time a highly significant census of environ- 
ments and structures up to redshift one, as well as a full morphological description of the galaxy population. In 
this paper we present a study aimed to constrain the evolution, in the redshift range 0.2 < z < L2, of the mass 
content of different morphological types and its dependence on the environmental density. We use a deep mul- 
ticolor catalog, covering an area of ^^0.7 D° inside the COSMOS field, with accurate photometric redshifts (i < 
26.5 and Az/iZspec+i) ~ 0.035). We estimate galaxy stellar masses by fitting the multi-color photometry to a 
grid of composite stellar population models. We quantitatively describe the galaxy morphology by fitting PSF 
convolved Sersic profiles to the galaxy surface brightness distributions down to F814 = 24 mag for a sample of 
41300 objects. 

We confirm an evolution of the morphological mix with redshift: the higher the redshift the more disk- 
dominated galaxies become important. We find that the morphological mix is a function of the local comoving 
density: the morphology density relation extends up to the highest redshift explored. 

The stellar mass function of disk-dominated galaxies is consistent with being constant with redshift. Con- 
versely, the stellar mass function of bulge-dominated systems shows a decline in normalization with redshift. 
Such different behaviors of late-types and early-types stellar mass functions naturally set the redshift evolution 
of the transition mass. 

We find a population of relatively massive, early-type galaxies, having high SSFR and blue colors which 
live preferentially in low density environments. The bulk of massive (> 7 x IO'^Mq) early-type galaxies have 
similar characteristic ages, colors, and SSFRs independently of the environment they belong to, with those 
hosting the oldest stars in the Universe preferentially belonging to the highest density regions. 

The whole catalog including morphological information and stellar mass estimates analysed in this work, is 
made publicly available. 

Subject headings: galaxies: evolution - galaxies: mass function - galaxies: statistics - galaxies: fundamental 
parameters - surveys 



> 

o 

00 
O 

o 



1. INTRODUCTION 

Galaxy formation and evolution has been a very actively 
debated topic of observational cosmology in the last years. 
Only recently models and observations are converging toward 
a unified and coherent picture. 

On one side the observations of the Universe, with some 
controversies on the way, have finally agreed that half of the 
present-day stars was already in place at z « 1 and that the 
same mass deficit affects, within the measured accuracies, 
the high-mass end: basically, the galaxy stellar mass func- 
tion evolves from the local determination to redshift one only 
by a normalization factor of about two. This means that, both 
at high and at low masses, there is roughly a factor of two 
difference between the comoving number density at redshift 

Based on observations made with the Advanced Camera for Surveys on 
board the NASA/ESA Hubble Space Telescope (GO Proposal 9822) 
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one and the local value (e.g. Dickinson et al. 2003, Rudnick 
et al. 2003, Drory et al. 2004a, Bundy et al. 2005, Pannella 
et al. 2006, Borch et al. 2006, Fontana et al. 2006, Pozzetti 
et al. 2007, Arnouts et al. 2007, Perez-Gonzalez et al. 2008, 
Cowie & Barger 2008, Drory & Alvarez 2008, Marchesini 
et al. 2008, Ilbert et al. 2009b). Measurements of the star for- 
mation rate density over cosmic times (e.g. Giavalisco et al. 
2004, Gabasch et al. 2004b, Bouwens et al. 2004, Sawicki & 
Thompson 2006, Bouwens et al. 2007) are also in good agree- 
ment with half of the present-day stars being already born at 
redshift L 

On the other side the models, linking the hierarchical 
growth of dark matter structures to the observed galaxy prop- 
erties by means of simplified prescriptions for the formation 
of baryonic systems, predict that galaxies form in a bottom-up 
fashion by following the cosmological destiny of dark mat- 
ter halos. Massive galaxies, in these models, assemble most 
(50%) of their stellar mass via merging at z < 1 (De Lu- 
cia et al. 2006), but most of these accretion events are red 
mergers (i.e. gas-free mergers that imply no induced star- 
formation), which are very difficult to detect from the obser- 
vational point of view (e.g. Hopkins et al. 2008, Cox et al. 
2008). The latest realizations of these models (e.g. Bower 
et al. 2006, Kitzbichler & White 2007, Cattaneo et al. 2008, 
Fontanot et al. 2009) have been able to fully reproduce the 
galaxy stellar mass function up to high redshifts, thus recon- 
ciling the theoretical bottom-up assembly of dark matter halos 
to the claimed top-down assembly of galaxies (e.g. Cimatti 
et al. 2006). 
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In hierarchical models star formation is supposed to take 
place only in disk structures and in gas-rich mergers (see also 
Dekel et al. 2009 for an up to date perspective). This latter 
event would exhaust the residual gas and create a dynami- 
cally hot system: a bulge galaxy. Also dynamical instabilities 
and minor mergers are able to destroy the disk structures and 
finally create an elliptical galaxy. Detailed predictions in the 
literature do not exist but the general expectation is that the 
fraction of bulge-dominated massive galaxies increases with 
time. 

The study of the evolution of the morphological mix in the 
galaxy stellar mass function is a decisive tool to put sensi- 
tive constraints on the models, because it could give unique 
insights in both the feedback and the merging processes in 
galaxy evolution (see for a detailed discussion Cole et al. 
2000). 

Galaxies have been often classified, both at low and at high 
redshift, as "red sequence" and "blue cloud" objects based on 
their broad band colors (Baldry et al. 2004, Fontana et al. 
2004, Giallongo et al. 2005, Bundy et al. 2006). However this 
kind of classification when used as a proxy for the galaxy mor- 
phology (see e.g. Mobasher et al. 2009), suffers from some 
obvious drawbacks: a disk galaxy populated by an old stellar 
population would be classified as an early-type object, and 
vice versa; a starburst galaxy, heavily obscured by dust, would 
be classified as a red, and dead, early-type object. 

Thanks to the advent of the Advanced Camera for Surveys 
(ACS; Ford et al. 2003) on board the Hubble Space Tele- 
scope (HST), it has become possible to investigate the optical 
rest-frame emission of galaxies up to redshift one and slightly 
beyond with sub-kiloparsec resolution and on relatively wide 
fields for hundreds of thousands of galaxies. 

Given the exponentially growing number of galaxies avail- 
able to be morphologically classified it has become urgent 
for galaxy evolution studies to build up robust and automated 
tools for morphological and structural classification that could 
reliably substitute the visual classification (see the discussion 
in Bamford et al. 2009 on unavoidable drawbacks affecting 
the automated approach). The efforts on the path to the auto- 
mated image classification have essentially split in two main 
branches, parametric and nonparametric methods, each hav- 
ing its own advantages and drawbacks. 

Parametric methods are based on the idea that galaxy light 
profiles may be sufficiently well described by analytic formu- 
las. Indeed parametric profiles have been shown to well de- 
scribe galaxy structural properties since many decades now. 
The law (de Vaucouleurs 1948) describing almost per- 
fectly luminous early-type galaxy hght profiles, and the ex- 
ponential law (Freeman 1970) which instead well fits the 
light profile of disk-dominated galaxies, have been exten- 
sively used to derive galaxy physical quantities, such as the 
effective radius and surface brightness. This has in fact al- 
lowed to formulate and study galaxy scaling relations, such 
as the TuUy-Fisher (TuUy & Fisher 1977) or the Fundamental 
Plane (Dressier et al. 1987, Djorgovski & Davis 1987, Bender 
et al. 1992), which have enormously boosted our knowledge 
and understanding of galaxy evolution and formation. 

The and the exponential laws may be used in the same 
modeling to describe both the bulge and disk photometric 
components of a single galaxy (see e.g. Simard et al. 1999). 
This approach is in fact doable only for very high signal to 
noise objects as the number of parameters involved becomes 
fairly large and hence the fitting itself heavily degenerate. For 



this reason it has become very popular the use of the more 
simple and robust single component Sersic (1968) model'. 
This is a convenient generalization of the R^ profile with 
where n is called the Sersic index. The variation of the index n 
allows to describe very different galaxy light profiles from the 
classical R^ (n = 4), to the late-type galaxy exponential profile 
(n = 1), to the massive elliptical galaxy profiles (n > 3.5) and 
down the unresolved point-like Gaussian profiles (n = 0.5). 

The Sersic profile fitting, thanks to its flexibility and robust- 
ness, has become a preferred tool in the last years to classify 
galaxies in broad morphological classes (e.g. Shen et al. 2003, 
Trujillo et al. 2004, Ravindranath et al. 2004, Barden et al. 
2005, Mcintosh et al. 2005, Pannella et al. 2006, Haussler 
et al. 2007), possibly combined to the "bumpiness" parameter 
as done in Blakeslee et al. (2006) and van der Wei (2008). 

The main idea behind the nonparametric approach is to 
measure model independent quantities, such as the concentra- 
tion, the asymmetry (Abraham et al. 1996, Conselice 2003), 
the M20 and the Gini parameters (Abraham et al. 2003, Lotz 
et al. 2004, Scarlata et al. 2007a) of the galaxy hght distri- 
bution and to correlate these quantities to the visual galaxy 
appearance. Albeit being model independent, and hence con- 
ceptually more effective in describing the variegated galaxy 
morphology zoo, these classification schemes are not free 
from severe biases as addressed in Blanton et al. (2003) and 
more recently in Lisker (2008). 

Having in mind to continue the study presented in this paper 
by investigating also the evolution of galaxy sizes, we opted 
for the parametric approach, that delivers both morphological 
classifications and structural parameters. 

In this context the COSMOS (Scoville et al. 2007b) sur- 
vey is the ultimate effort to cover a sensitive patch of the sky 
(2a°) with HST, going factors of 10/100/1000 wider than pre- 
vious efforts like the HDF-N (Williams et al. 1996), GEMS 
(Rix et al. 2004), the Great Observatories Origins Deep Sur- 
vey (GOODS; Giavalisco et al. 2004), the FORS Deep Field 
(FDF; Heidt et al. 2003) among many others. 

The high resolution imaging has been complemented in 
COSMOS by a state-of-the-art multiwavelength coverage all 
the way from the X-ray to the radio. By observing the Uni- 
verse in all its colors COSMOS is becoming a real break- 
through in observational cosmology. The need to tackle and 
understand the conspiracies of Nature embedded in the Uni- 
verse variance makes of this wide and deep survey a unique 
and unprecedented tool for galaxy evolution studies. 

In this work we rely on a deep and accurate (i < 
26.5 and Az/{zspec + ^)^ 0.035) photometric redshift cata- 
log (Gabasch et al. 2008) complemented with a morpholog- 
ically classified sample (F814 < 24) to study the contribution 
of galaxies of different morphologies to the redshift evolution 
of the stellar mass density, as well as its dependence on the en- 
vironmental density, over an area of ~0.7 in the COSMOS 
field. 

Brinchmann & Ellis (2000), Bundy et al. (2005), Pannella 
et al. (2006) and Franceschini et al. (2006) have carried out 
studies on the evolution of the galaxy stellar mass function 
split by morphological types, in a way qualitatively similar to 
the one we present here, but based either on shallower samples 
and visual morphological classification or on much smaller 
areas heavily affected by cosmic variance. Likely due to the 
small surveyed area, no one of the just mentioned studies de- 

' We refer to Appendix A for the analytic expression of the S&sic profile. 
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Fig. 1. — Redshift distribution of tlie 41300 galaxies belonging to the 
morphological catalog used in the present work. Overdensities at redshifts 
0.4, 0.7 and 1 are quite prominent, as well as the underdensities at redshifts 
0.5 and 0.8. 

fined environmental properties of galaxies in order to contrast 
the effect of local environment against cosmic time on mor- 
phological evolution. 

This paper is organized as follows: in § 2 we discuss the 
ground-based dataset and the photometric redshifts on which 
this work is based, in § 3 we discuss how galaxy stellar masses 
have been estimated, in § 4 we present the ACS data used in 
this work, the source extraction and the galaxy number counts, 
in § 5 we briefly describe the quantitative morphological anal- 
ysis, referring to a more detailed discussion in Appendix A, in 
§ 6 we present the evolution of the morphological mass func- 
tion and mass density, in § 7 we discuss the effect of local 
environments on galaxy morphological evolution and finally 
in § 8 we summarize our results and draw our conclusions. 

Throughout this work, we use AB magnitudes and 
adopt a A cosmology with ^Im = 0.3, J7a = 0.7, and 
//o = 70kms-iMpc"'. 

2. GROUND-BASED DATA AND PHOTOMETRIC REDSHIFTS 

The ground-based data used in this paper combine the uB- 
VrizK COSMOS dataset ' with proprietary imaging in the H 
band. A highly homogeneous multi-wavelength i-band se- 
lected catalog is used to derive accurate photometric redshifts. 
A full description of the NIR data acquisition, data reduction 
and catalog assembling, and of the photometric redshifts es- 
timation, is presented in Gabasch et al. (2008). Here in the 
following, for the sake of clarity, we give only a brief sum- 
mary. 

Our cleaned i-selected catalog- comprises 293 377 objects 
over an area of '--^0.7 square degrees down to a limiting mag- 
nitude ^-^26. 5. There is a good agreement between literature 
data and this catalog number counts up to the 25th magni- 
tude. At the bright end (i<21) the literature number counts 
are sensibly higher since most of the very bright objects are 

' publicly available at http://irsa.ipac.caltech.edu/data/COSMOS 
^ publicly available at http://www.mpe.mpg.de/~gabasch/COSMOS 
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Fig. 2. — The distribution of galaxy stellar masses, derived using a Salpeter 
(1955) initial mass function, as a function of redshift for the morphological 
catalog used in this work (empty squares). The sohd blue line shows the 
estimated mass completeness of the morphological sample at each redshift 
(see text for details). For comparison, data from the morphological catalog 
used in Pannella et al. (2006) are overplotted (yellow and magenta symbols 
for the FDF and GOODS-S fields, respectively). 



saturated in the Subaru i-band images and thus not present in 
the final catalog. Since objects brighter than i=21 are likely 
to be important for the determination of the high-mass end of 
the galaxy stellar mass function at low and intermediate red- 
shifts, we added those objects (1896) back to our catalog by 
complementing it with the COSMOS legacy photometric cat- 
alog described in Capak et al. (2007b), where CFHT i band 
photometry was used for this bright tail. 

Photometric redshifts were derived using the technique de- 
scribed in Bender et al. (2001), Gabasch et al. (2004a) and 
Brimioulle et al. (2008). In short, the method consists of : 
i) checking photometric zeropoints and, if necessary, deter- 
mining photometric offsets by comparing theoretical and ob- 
served stellar locii, ii) computing object fluxes in a fixed aper- 
ture (2.0") from seeing-matched images, iii) determining a 
redshift probability function P(z) for each object by matching 
the object's fluxes against a set of template spectra covering a 
wide range of stellar population ages and star-formation his- 
tories. 

The additional use of the GALEX FUV and NUV bands 
(Zamojski et al. 2007) allowed some photo-z degeneracies to 
be broken, and a final accuracy of /S.z/{Zspec+ 1) ~ 0.035 to 
be reached. Figure 1 shows the redshift distribution of the 
subsample with available morphological information used in 
this work.^ 

The morphological catalog, including stel- 
lar mass estimates, analysed in this work is 
made publicly available at the following address: 
http ://w WW. aoc. nrao . edu/~^mpannell/data. html . 

^ In the last stages of this work the COSMOS Legacy released a new cat- 
alog of spectroscopic redshifts within the zCOSMOS DR2 (Lilly et al. 2007) 
as well as an updated catalog of photometric redshifts (Ilbert et al. 2009a). 
We refer to Appendix B for a comparison with these two catalogs. 



4 



Pannella et al. 




TABLE 1 

Log Mt/M© COMPLETENESS VERSUS REDSHIFT 
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Fig. 3. — Cosmic stellar mass density up to redshift 1 for objects with 
Log M, >^ 10.8, corresponding to the mass completeness value for the highest 
redshift considered. Horizontal dashed lines show the local mass density 
from Bell et al. (2003) after correcting for the different IMF and applying 
the same mass cut. The error bars represent the cosmic variance contribution 
(Somerville et al. 2004). Assuming that the mass function evolves only in 
normalization over this redshift range, adding 0.28 dex to the plotted values 
gives the total mass densities. 



3. COMPUTING MASS-TO-LIGHT RATIOS 

The method we use to infer stellar masses from multi- 
wavelength photometry is described in detail in Drory et al. 
(2004b). It is based on the comparison of the observed multi- 
color photometry with a grid of stellar population synthesis 
models produced with the Bruzual & Chariot (2003) code. 

We parametrize the possible star formation histories (SFHs) 
by a two-components model, consisting of a main, smooth 
component described by an exponentially declining star for- 
mation rate tp(t) oc exp(-j/r), linearly combined with a 
secondary burst of star formation. The main component 
timescale t varies in £ [0.1, oo] Gyr, and its metallicity is 
fixed to the solar metallicity. The age of the main compo- 
nent, t, is allowed to vary between 0.5 Gyr and the age of the 
Universe at the object's redshift. 

The secondary burst of star formation is modeled as a 
100 Myr old constant star formation rate episode of solar 
metallicity. 

We adopt a Salpeter (1955) initial mass function for both 
components, with lower and upper mass cutoffs of 0.1 and 
100 M0. 

Additionally, adopting the Calzetti et al. (1994) extinction, 
both the main component and the burst are allowed to ex- 

1 2 

hibit a variable amount of attenuation by dust with Ay S 
[0, 1 .5] and [0, 2] for the main component and the burst, re- 
spectively. This takes into account the fact that young stars 
are found in dusty environments and that the star light from 
the galaxy as a whole may be reddened by a (geometry de- 
pendent) different amount. We compare this set of mod- 
els with multi-color photometry of each object, computing 
the full likelihood distribution in the 5-dimensional parame- 
ter space (t,? ,Ay,/3,A^), the likelihood of each model being 
oc exp(-x^/2). 

To compute the likelihood distribution of mass-to-light ra- 
tios M/L, we weight the M/L of each model by its hkehhood 
and marginalize over all parameters. The uncertainty in M/L 
is obtained from the width of this distribution. 

For young objects, with relatively high burst fractions, the 
width of the M/L distribution is usually much wider in the op- 



tical than in NIR bands, while for quiescent objects the width 
of the M/L distribution is very similar in all bands. On av- 
erage, the width of the likelihood distribution of M/L at 68% 
confidence level is between ±0.1 and ±0.2 dex (using the B 
band M/L). The uncertainty in mass has a weak dependence 
on mass itself (increasing with lower S/N photometry), and it 
mostly depends on the spectral type: quiescent galaxies have 
more tightly constrained masses than star-forming ones. 

We estimated the galaxy stellar masses for all the objects in 
the i band selected catalog down to the limiting magnitude. 

In Figure 2 we show the distribution of galaxy stellar 
masses as a function of redshift for the morphological catalog 
used in this work. For comparison, and to show the impres- 
sive statistic available, we overplot the morphological cata- 
logs used in Pannella et al. (2006) from the FDF (yellow) and 
GOODS-S (magenta) fields. 

In Figure 3 we show the stellar mass density up to redshift 
1 for objects with Log > 10.8, which corresponds to the 
mass completeness value for the highest redshift considered. 
Dashed lines represent Itr estimates of the local mass density 
from Bell et al. (2003) after applying the same mass cut. The 
error bars represent the cosmic variance contribution, as esti- 
mated in Somerville et al. (2004), that is a factor 10 larger than 
statistical errors. Assuming that the mass function evolves 
only in normalization over this redshift range, one should add 
0.28 dex to the plotted values to recover total mass densities. 

The mass completeness at different redshifts was estimated 
as the mass of a maximally old stellar population"* having, at 
the redshift considered, an observed magnitude equal to the 
catalog completeness magnitude (e.g. Dickinson et al. 2003). 
In table 1 we provide the mass completeness values as a func- 
tion of redshift for the morphological catalog (F814W <24), 
as well as for the global photometric catalog (i < 26). 

4. HST ADVANCED CAMERA FOR SURVEYS IMAGING 

The 2U° COSMOS field was imaged in the F814W ACS 
filter for a total of 581 orbits. Within each orbit, four equal 
length exposures of 507 sec duration each (2028 sec total) 
were obtained in a 4 position dither pattern, designed to shift 
bad pixels and to fill in the 90 pixel gap between the two 
ACS CCD arrays. Adjacent pointings in the mosaic were po- 
sitioned with approximately 4% overlap in order to provide 
at least 3 exposure coverage at the edge of each pointing and 
4 exposure coverage over approximately 95% of the survey 
area. This multiple exposure coverage with ACS provided ex- 
cellent cosmic ray rejection. 

A full description of the ACS data processing including 
drizzling, flux calibration, registration and mosaicing is pro- 
vided in Koekemoer et al. (2007) and in Scoville et al. 
(2007a). 

The ACS images released to the public^ are sampled with 
0.05" pixels. The measured FWHM of the PSF in the ACS 
F814W band filter is w 0.11". 

^ We used a dust-free, passively evolving stellar population model, ignited 
by an instantaneous burst of sub solar (Z=0.008) metaUicity at z = 10. 
' Publicly available at http://irsa.ipac.caltech.edu/data/COSMOS/ 
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TABLE 2 

Galaxy number counts not corrected for incompleteness in 

THE FS14W PASSBAND (~0.7 D"). logiV AND cr\agN ARE GIVEN, WHERE 

iV IS IN UNITS OF 0.5mag"' deg"^ . 





FSUW 




mag 


log/V 




17.01 


1.576 


0.0969 


17.51 


1.801 


0.0728 


18.01 


2.137 


0.0481 


18.51 


2.416 


0.0343 


19.01 


2.655 


0.0258 


19.51 


2.867 


0.0201 


20.01 


3.047 


0.0163 


20.51 


3.252 


0.0128 


21.01 


3.429 


0.0104 


21.51 


3.63 


0.0082 


22.01 


3.793 


0.0068 


22.51 


3.974 


0.0055 


23.01 


4.148 


0.0045 


23.51 


4.318 


0.0037 


24.01 


4.501 


0.0030 


24.51 


4.671 


0.0025 


25.01 


4.802 


0.0021 


25.51 


4.872 


0.0020 


26.01 


4.769 


0.0022 


26.51 


4.01 


0.0053 




Fig. 4. — Galaxy number counts for the ACS/F814W passband are plotted 

as solid dots. Error bars account for poissonian errors only, and no correction 
for incompleteness has been applied. Literature data are plotted with other 
symbols and are taken from: Pannella et al. (2006) (asterisks, FDF), Leau- 
thaud et al. (2007) (open squares, COSMOS) and Benftez et al. (2004) (diag- 
onal crosses, VV 29). 



4.1. Source extraction, star-galaxy classification and 

cataloging 

Source extraction was performed with the SExtractor code 
(Berlin & Arnouts 1996) on each ACS tile in a cold manner, in 
order to minimize the artificial source splitting and maximize 
the number of pixels assigned to each object. 

We used 348 ACS tiles, namely only those overlapping with 
the ground based photometry used in this work. 

Taking advantage of the ACS high resolution imaging we 
could remove the point-like (stellar) objects using the struc- 
tural parameters output by SExtractor. We used the full width 
half maximum (FWHM), the half Ught radius, the neural net- 
work stellarity index and the total magnitude for each object 
to select point-like sources. 

In Figure 4 we show F814W band galaxy number counts, 
not corrected for incompleteness, compared with literature 
counts (Pannella et al. 2006, Leauthaud et al. 2007, Bemtez 
et al. 2004), and in table 2 we provide the plotted values. The 
ACS galaxy catalog has been cross-correlated with the ground 
based catalog with a 1 arcsec matching radius. We removed 
from the final catalog all the objects, a few percent of the total, 
for which one ground based entry had more than one entry in 
the ACS catalog. The final morphological catalog used in the 
following contains 41300 objects over ^0.7 and results to 
be more than 90 % complete down to F814W = 24. 

5. THE MORPHOLOGICAL ANALYSIS 

We use the package GIM2D (Simard et al. 1999) to fit 
PSF convolved Sersic (1968) profiles to the two-dimensional 
surface brightness distribution of each object, for all sources 
down to a limit of F814W=24. The PSFs used to convolve 
the profiles were obtained for each individual tile by median 
stacking about 50 high S/N stars. 

We tested our results by running also the GALFIT (Peng 
et al. 2002) code on some of the COSMOS ACS patches. The 
results from the two different codes are in excellent agree- 



ment, thus confirming both the robustness of our modelUng 
and the choice of flux limits. 

We further tested our results by running the G1M2D code on 
simulated images. This test was performed by adding to the 
real image, in a blank, pure-sky region, fake objects with Ser- 
sic profiles and structural parameters spanning a wide range 
of plausible values. The whole procedure is then run exactly 
as for real objects: we find that the profile parameters of sim- 
ulated objects are well recovered down to F814W=24. High 
Sersic index objects have very extended wings that, depend- 
ing on the total flux, can fall under the sky surface bright- 
ness, thus for these objects a lower Sersic index, total flux and 
effective radii are usually recovered (see also Sargent et al. 
2007). We point out that the same trend was found in Tru- 
jillo et al. (2006b) with a different fitting code, thus confirm- 
ing the "physical bias" as contrasted to a code failure. In the 
local Universe a tight correlation between Sersic index and 
luminosity/mass of giant elliptical galaxies (see e.g. Caon 
et al. 1993 and Graham et al. 2001) has been robustly estab- 
hshed. Due to flux limits and mass completeness, morpho- 
logical analysis of high redshift ellipticals is mostly limited 
to the most massive galaxies which one would expect to have 
relatively high Sersic index. Our simulations show that for 
such high Sersic index objects at faint apparent magnitudes 
(low S/N) we would recover Sersic indexes and effective radii 
biased toward lower values. While we will deal in detail with 
the galaxy size evolution in a future work, we argue here that 
the strong size evolution of early-type galaxies with redshift, 
recently reported in the literature (e.g. Trujillo et al. 2006a, 
Cimatti et al. 2008, van Dokkum et al. 2008), could be partly 
due to such bias. We try to cope with this "biased" Sersic in- 
dex by adopting the procedure described in the following. We 
perform extensive Monte-Carlo simulations to take into ac- 
count the Sersic index uncertainty (Ri30%, estimated through 
simulations) on our results. One thousand simulations of 
the morphological catalog were generated. Objects with Ser- 



6 



Pannella et al. 



sic index smaller than 2.5 were perturbed within a Gaussian 
of sigma equal to the 30% of the measured value. Objects 
with Sersic index larger than 2.5 were perturbed within only 
the positive side of the Gaussian so that their value could only 
get larger than the measured ones. We tried also to implement 
different solutions to this correction procedure and the final 
results were always in good agreement. More details on the 
morphological analysis are given in the Appendix A. 

6. THE EVOLUTION OF THE GALAXY STELLAR MASS FUNCTION 
BY DIFFERENT MORPHOLOGIES 

Pannella et al. (2006) show that there is a good correlation 
between the average visual and automated classifications as 
parametrized by the morphological type T, assigned accord- 
ing to the de Vaucouleurs et al. (1991) classification scheme, 
and the Sersic index n^er, respectively (see Figure 1 in that pa- 
per). Following Pannella et al. (2006), we split our sample, 
according to the Sersic index value in three broad morpholog- 
ical classes: early-type (n > 3.5, which corresponds to bulge- 
dominated/elliptical galaxies), intermediate (2 < n < 3.5) and 
late-type (n < 2, which corresponds to disk-dominated galax- 
ies). 

We restrict our analysis to z < 1.2 to limit the bias intro- 
duced by the rest-frame emission band shifting, namely the 
effect of the morphological k-correction. A number of stud- 
ies in the last years have shown, although not in a conclusive 
and quantitative way, that this is indeed a fair assumption (e.g. 
Scarlata et al. 2007b, Abraham et al. 2007). 

We performed extensive Monte-Carlo simulations to take 
into account the effect of mass uncertainties (« 0.2dex) on our 
results. One thousand simulations of the mass catalog were 
generated, perturbing each mass within a Gaussian of sigma 
equal to its error. Unless stated differently in the relevant Fig- 
ure captions, we use the median values of the Monte-Carlo 
simulations in all Figures. Error bars take into account both 
poissonian errors on the median counts (Gehrels 1986), and 
16-84''' percentile values of each distribution. 

In Figure 5 we show the VNmax corrected number density 
evolution split by morphological type, up to redshift 1.2, for 
objects in stellar mass ranges Log M*=[10.6-10.8], [10. 8-11], 
[11-11.2], and [11.2-12], as indicated in the upper right cor- 
ner of each panel. By multiplying the number densities by 
the average mass in the corresponding mass range, each panel 
can be interpreted as the evolution of the total mass densities 
contributed by objects in that stellar mass range. The dotted 
vertical lines show the estimated redshift completeness for the 
relevant mass bin. 

We find that early-type objects always dominate the high- 
mass end. The lower the stellar mass, the more late-type 
objects tend to dominate the relative contribution to the to- 
tal number/mass density. The same trend is true by looking 
at increasing redshift: the higher the redshift the more the rel- 
ative contribution of early and late types gets closer. For the 
lowest mass range [10.6-10.8] the transition from a bulge- 
dominated to a disk-dominated stellar mass budget happens 
around redshift 0.7; beyond that redshift late-types are always 
contributing more to the total mass budget. In the mass range 
Log M» = [10.8-1 1] the transition is occurring at z « 1.1. 

At each redshift, a transition mass can be identified where 
the transition from a bulge-dominated to a disk-dominated 
stellar mass budget takes place. Based on the very good statis- 
tics available in this work, we are able to confirm the main 
conclusions of Pannella et al. (2006): i) the morphologi- 
cal mix at the high-mass end evolves with redshift; ii) the 
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stellar mass bin is complete. Symbols are colored according to the morpho- 
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transition mass increases with redshift. At z ~ 0.7 we find 
a value of the transition mass approximately consistent with 
the local value M « 5 x 1O'"M0, as measured in Bell et al. 
(2003) by using the concentration parameter to discriminate 
between early and late type objects in a complete sample ex- 
tracted from local surveys. At z 1 the late and early types 
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contributions become comparable at M « 1 x 1O"M0. 

In a fixed mass range the number densities of late-types are 
consistent with being constant with redshift, while the early 
types follow the more general total declining. At a fixed red- 
shift, the number density of early types of different masses 
is very similar (in +- 0.1/0.2 dex), while the late type num- 
ber densities move by one order of magnitude passing from 
the high to the low stellar mass range. This naturally sets the 
evolution of the transition mass with redshift. 

The slope of the redshift evolution of total number densities 
is consistent in all mass ranges, at least up to the common 
redshift completeness, meaning that the mass function in its 
high-mass end, i.e. at least down to M « 4 x lO'^^M©, does 
not evolve in shape but only in normalization with redshift. 

7. DISENTANGLING THE ROLE OF ENVIRONMENT ON 
GALAXY MORPHOLOGICAL EVOLUTION 

In Figure 5 there is obvious scatter around the mean number 
density decUne in all the mass bins. However, an underden- 
sity at z^O.5 shows up quite strikingly. The deficit of objects 
is also quite evident in the global redshift distribution of our 
morphological catalog, as shown in Figure 1. It gets more 
and more significant with increasing stellar mass. The deficit 
of galaxies at this redshift is a factor of three for the highest 
mass objects, and less than a factor of two for objects in the 
lowest mass bin explored. This differential deficit is in qual- 
itative agreement with galaxy evolution dependence on the 
local density as expected in current galaxy formation models 
(see Lee 2006 and references therein): the more underdense 
is the environment the later massive objects will assemble be- 
cause there is in general a time delay in starting star formation 
and also there are less small galaxy units to assemble the mas- 
sive giant galaxies. 

From the morphological point of view, the evolution in this 
underdense redshift bin is delayed at all stellar masses as com- 
pared with the contiguous redshift bins. The early-type frac- 
tion with respect to both the total number density and the late- 
type number density is much lower than those in the contigu- 
ous redshift bins, in agreement with the expectations from the 
morphology-density relation. We also note that in the high- 
est mass bin, in this underdense region there are no massive 
late-type galaxies. 

We notice that the morphological mix in this underdense 
redshift bin is, at all masses probed, very close to the morpho- 
logical mix at zwl.l, which has a comparable total number 
density of objects. 

Taking advantage of the good statistics available in this 
study, we wish to further check the importance and role of 
environment on galaxy morphological evolution. 

Estimating the local volume density in a photometric red- 
shift survey is not straightforward, and we try here to make it 
as simple and unbiased as possible. 

First of all one has to decide which are the objects tracing 
the overdensities in a way not biased by the flux limited sur- 
vey: a density-defining population (DDP, Croton et al. 2005) 
of galaxies is needed. The objects have to be bright enough 
to be seen in the whole redshift range considered. Dealing 
with masses and morphologies, one should also take care of 
the selection effects introduced by the different M/L ratios 
associated with different steUar populations. BasicaUy, old 
stellar populations, which more often but not always reside in 
bulges, have much higher M/L ratios with respect to young 
stellar populations, often but not always populating disks. A 
pure luminosity cut - even taking into account passive evolu- 
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Fig. 6. — Number density evolution for objects with Log M, > 10.8. 
The solid line, p(z) = -2,62 — 0,36 X z, shows the Unear fit to the filled dots 
(i,e, the number densities obtained over the whole field in redshift bins of 
Az = 0,2), The dashed line sets the l.4xp(z) limit above which we select 
the overdense regions. Below the dashed-dot (0,5 X p(z)) line we select the 
underdense regions, while in between the two lines we select the medium- 
dense regions. The empty squares overplotted show the number densities in 
all the cells we have split our survey volume. 



tion - would thus preferentially cut out of the sample high 

M/L objects (early-types) with respect to low M/L objects 
(late-types). To avoid this kind of bias, we use as DDP all 
objects with masses down to the mass completeness value at 
the highest redshift we are exploring (7x10^'^ solar masses). 

Another issue is the relative scarcity of these objects (being 
in the exponential cutoff of the mass function), so one has to 
avoid reaching the shot noise regime when calculating number 
densities. 

With this in mind, we proceed as follows: we first fit the 
number density evolution with redshift for objects more mas- 
sive than 7x10"* solar masses, obtaining in this way an ex- 
pected value of the number density, p{z), at each redshift z. 
We then split the entire redshift range in slices of Az = 0.1 
(corresponding to differential depths of about about 327 co- 
moving Mpc at redshift 0.5 and 244 at redshift 1). This slicing 
along the redshift axis turns out to be quite robust against the 
photo-z errors, in fact this corresponds to about 2 times the 
formal photo-z error at redshift 0.5 and 1 .5 times at redshift 1 . 
We then split the entire survey volume in cells of 2.8 comov- 
ing Mpc by side at each redshift slice, that is the maximum 
contiguous angular area we have available in the lowest red- 
shift bin. Finally, we calculate the comoving number density 
of each cell at the different redshifts and assign the cell galaxy 
population to one of three density classes: underdense (piz) < 
0.5 X p(z)), medium-dense (0.5 x p(z) < p(z) < 1 .4 x p(z)) and 
overdense(p(z) > 1.4 x piz)). The Umits used here to define 
environmental classes are chosen to balance the needs to map 
sensibly different environments and to assign a relatively large 
number of objects to each of the three classes. 

In Figure 6 we show the number density evolution for ob- 
jects with Log M, > 10.8. The sohd Une, p(z) = -2.62- 
0.36 X z, shows the linear fit to the filled dots (i.e. the number 
densities obtained over the whole field in redshift bins of Az 
= 0.2). The dashed line sets the 1.4xp(z) limit above which 
we select the overdense regions. Below the dashed-dot line 
(0.5 X p{z)) we select the underdense regions, while in be- 
tween the two lines we select the medium-dense regions. The 
empty squares overplotted show Ihc number densities in all 
the cells we have spht our survey volume. 

We run 100 realizations of the procedure adopted to define 
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Fig. 7. — Galaxy stellar mass functions as a function of redshift, morphology and local environment. Black symbols refer to total values. The vertical axis is 
in units of [h^^Mpc^^^dex"']. A vertical dotted line indicates the mass completeness limit in each redshift bin. Colored symbols refer to morphological classes, 
see Figure 5 for details. Dashed colored lines are only intended to guide the eye. Solid black lines show the local filed mass function determinations from Cole 
et al. (2001) and Bell et al. (2003). 



environmental densities by randomizing the initial position of 
the grid in the sky and also by adopting different solutions for 
the redshift slicing of the total volume. The outcome of this 
testing convinced us that the results we are going to show and 
discuss in the following are robust and not biased against the 
specific procedure. 

7.1. Morphological evolution and mass growth of 
early and late-type galaxies in 0.25 < z < 1.15 

In Figure 7 we show the galaxy stellar mass function as a 
function of redshift, split by morphology and for the differ- 
ent environmental densities as defined in the previous section. 
Each of the panels contains about 2000 objects. Vertical dot- 
ted lines set the estimated mass completeness at the different 
redshifts. For reference in all panels are overplotted local es- 
timates for the stellar mass function (Cole et al. 2001, Bell 
et al. 2003), which are fixed and independent from the envi- 
ronment. They are only meant as a reference to look at differ- 
ential evolution in different environments. 



In the lowest range explored, at z ^ 0.4, the mass com- 
pleteness reaches masses as low as 7x10'^Mq. This enables 
us to explore the behavior of the slope of the mass function 
in different environments. The slope of the mass function is 
marginally evolving from the lower to the higher density re- 
gions, becoming increasingly shallower and consistent with 
the local field one. This trend is in fact in good agreement 
with what found in the local Universe (Baldry et al. 2006). 

The transition mass moves from Log ~ 10.6 at high 
density to Log M,^ ~ 10.9 at low density. 

In the overdense and medium-dense environments the 
shape of the mass function is in a very good agreement with 
the local determinations. In the overdense regions, this hap- 
pening at all redshifts, a very massive bulge-dominated popu- 
lation seem to appear that overcomes the local field determi- 
nation. 

Figure 7 shows that the morphologically split mass func- 
tion evolves along both the redshift and the density axes. The 
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evolution is differential though: the redshift evolution is much 
milder than the environmental one. 

This is in fact a crucial point in interpreting the results since 
in almost all previous studies there has been no chance to 
split these two very different contributions. This is better ex- 
plained by the following examples: the total mass function, 
and the morphological mixtures are basically the same, up 
to the common mass completeness value, for overdense vol- 
umes at redshift 1 and medium-dense volumes at redshift 0.4. 
The transition mass is close to Log M, = 11 at redshift 1 in 
medium-dense volumes and it is already at that mass in the 
underdense environments at z ~ 0.4. This shows how by an- 
alyzing small patches of the sky one could mix up different 
environments at different redshifts and hence get the wrong 
evolutionary paths for different morphologies. 

With the aim of trying to better quantify the morphological 
mix evolution as a function of redshift and environment, we 
integrate the mass function over the common complete mass 
range in the three redshift bins and in the three different envi- 
ronments. 

In Figure 8 we show both the early-type fraction (lower 
panel) and the early-to-late type classification ratio (upper 
panel) at different redshifts and for different environments 
(low-cyan, medium-magenta, black-high density). 

At all redshifts a morphology-density relation, as already 
reported in the literature (e.g. Smith et al. 2005), is showing 
up: denser volumes contain also a higher relative fraction of 
bulge-dominated galaxies with respect to lower density vol- 
umes. Although it is difficult to compare with previous re- 
sults given the fact that in our study for the first time, to our 
knowledge, a mass complete sample has been adopted, it is 
very reassuring that qualitatively we are in agreement with 
previous studies. This confirms the robustness of our galaxy 
environment definition procedure. 

As a function of redshift, in the medium and low density 
environments there is an increase of the early-type population 
with cosmic time which is more pronounced between redshift 
1 and 0.7 and then flattens out to redshift 0.4. In the high 
density environments the early-type fraction grows between 
redshift 1 and 0.7, and then seems to turn-over down to red- 
shift 0.4. A similar behavior has been actually found in Capak 
et al. (2007a). In our opinion the sky patch, we are using here, 
might not cover enough cosmic volume at redshift around 0.4 
to allow an exact environmental definition. This shows up in 



Figure 6 where all but two sub-volumes are in fact assigned 
to the medium dense environment bin. 

7.2. The ages of early-types in dijferent environments 

In Figure 7.2 we show the ages of massive early-type galax- 
ies as a function of redshift and for different environments 
as obtained from the SED fitting procedure outlined above. 
The redshift ranges considered are indicated in the top label 
of each panel. We cut the total sample to the common com- 
pleteness limit (> 7x 10'"Mq). The dotted, dashed and solid 
histograms represent the high, medium and low density envi- 
ronments, respectively. In each panel, the difference in bil- 
lion years between the mean ages in the two most extreme 
environments (the high and low density regions) is indicated. 
Such differences deviate from zero at a two sigma level at 
all redshifts, and might suggest an earlier formation epoch of 
early-type galaxies in high density environments with respect 
to those living in a low density region. The results found are in 
good agreement with recent studies (e.g. Thomas et al. 2005, 
van Dokkum & van der Marel 2007, Gobat et al. 2008, Rettura 
et al. 2008). 

From Figure 7.2 three things are worth noting: first, the 
histograms are quite similar to each other within each single 
panel (as confirmed by a Kolmogorov-Smirnov test) but for 
the oldest and youngest tails of the distribution which prefer- 
entially belong to the highest and the lowest density regions; 
second, the width of the distributions is increasing with the 
cosmic clock, this being partially driven by the narrower range 
of ages allowed at earlier cosmic times; third, in the lowest 
redshift bin a number of, relatively young, early-type objects 
are present which cannot be the descendants of the higher red- 
shifts early-type populations simply evolved by passive evo- 
lution. 

We interpret these evidences as follows: i) young massive 
early-type galaxies are preferentially found in low density en- 
vironments at all redshifts, ii) the bulk of the early-type stel- 
lar populations have similar characteristic ages independent 
of the environment; iii) the oldest objects in the Universe, or 
better saying the oldest stars in the Universe, happen to belong 
to its highest density regions; iv) newly formed, or somehow 
rejuvenated, bulge-dominated galaxies are entering the sam- 
ple with cosmic time. This latter finding might explain why 
the difference in age seems to increase with the cosmic clock 
and reach the 2 Gyrs quoted by Thomas et al. (2005) in the 
local Universe. 

The estimate of stellar population ages suffers from some 
important degeneracies and has to be treated with a pinch 
of salt. The SED modelling we have described in sec- 
tion 3 only considers stellar population models with solar 
metallicity. This means we are not tackling the well known 
age-metallicity degeneracy: old stellar populations with low 
metallicities and younger stellar populations with high metal- 
licities may well have very similar SEDs. Only recently 
Cooper et al. (2008), but see also Mouhcine et al. (2007) 
for different conclusions, reported on a mild observed cor- 
relation between metallicity and environment: massive galax- 
ies in high-density environments are slightly more metal rich 
than galaxies belonging to low-density environments. If such 
a correlation is indeed at work our results are biased and the 
correct differences in age would be even smaller than the ones 
reported in the panels of Figure 9. 

We remind that our classification scheme does not al- 
low any separation between ellipticals and early-type spirals. 
Therefore, as happens in the local Universe (see e.g. Dressier 
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Fig. 9. — Ages of massive (> 7x IO'^Mq) early-type objects at different redshift and for different environments. The difference, in billion years, between the 
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1980) our early-type samples will be likely dominated by el- 
liptical and SO galaxies in overdense environments and by 
early-type spirals in low density regions. 

7.3. Red sequence and blue cloud 

In Figure 10 we show the rest-frame U-B color as a function 
of mass, redshift, morphology and environment. Morphologi- 
cal classification is highlighted in different colors (red = early, 
green = intermediate, blue = late). 

Contours are used to better show the different behavior of 
the three morphological classes. The contours correspond 
to equally spaced levels between the minimum and maxi- 
mum density of the data points of each morphological set. 
The first three rows of panels starting from the top, show the 
high, medium and low density environments. The bottom row 
shows instead the total population, i.e. all the different en- 
vironments together. Redshift ranges are indicated, for each 
column, in the top label. 

The red-sequence, mostly populated by early-types, in all 
environments becomes slightly redder with cosmic time be- 
cause of the ageing of the stellar populations. The evolution 
in the rest-frame U-B color is mild, as also shown for instance 
by Willmer et al. (2006) up to zwl.3. The red-sequence up- 
per envelope reddens by A(U-B) x 0.25 mag in the redshift 
range 0.25< z <1.15, in agreement with simple passive evo- 
lution expectations. 

The bulge-dominated galaxy population is mostly confined 
to the red-sequence with some peculiarities dependent on en- 
vironments. At all redshifts, in higher density environments 
early-types tend to be more confined to the massive end of the 
red- sequence. On the other hand, in lower density environ- 
ments early-types tend to populate the red-sequence down to 
lower masses (as it is easier to see at low redshifts) and reach 
into the blue-cloud. This definitely points toward a difference 
in age between the bulge-dominated populations in different 
environments. 

It is worth noting the evolution of the intermediate-type 
population: these are mainly spiral galaxies hosting fairly 
massive central bulges. At redshift 1, intermediate-type 
galaxies overlap in the color-mass diagram with both the 



early and late populations in all environments. As time goes 
by, this class of objects tend to leave the very high mass end 
of the red sequence in favour of bulge-dominated galaxies, 
and to drift toward the bluer late-type population. This sug- 
gests that part of these intermediate-type galaxies have moved 
to the early-type population. The change of this class of 
objects with cosmic time correlates with the environmental 
density: at low redshifts, intermediate-type objects in more 
dense environments tend to occupy the valley between the 
red-sequence and the blue-cloud, while in low density en- 
vironments they split in two populations, with the majority 
well overlapping the late-type population, and a minor sub- 
population remaining in red-sequence together with the mas- 
sive bulge-dominated galaxies. 

The late-type population is confined at all redshifts and in 
all environments to the blue cloud. Given the mass complet- 
ess limits it is definitely hard to make strong conclusions, but 
these data suggest the presence, possibly more pronounced at 
higher redshift and in lower density environments, of a pop- 
ulation of massive disk-dominated galaxies reaching toward 
the red-sequence. 

In section 7.1 we have investigated the evolution with red- 
shift of the morphology-density relation (see Figure 8). Since 
Figure 10 shows that at all redshifts early-types populate 
mainly the red sequence while late-type galaxies are mostly 
in the blue cloud, we can infer from our data that a color- 
density relation is in place up to the highest redshift explored. 

Cucciati et al. (2006) analysed a sample of 2900 galaxies 
with secure spectroscopic redshift determinations drawn from 
the VIMOS-VLT Deep Survey (VVDS, Le Fevre et al. 2005) 
in an area 4 times smaller than the one considered in this 
work and with «32% sampling, slightly biased toward faint 
objects, of the whole galaxy population. They obtain a very 
similar result to the one we present in this work up to redshift 
1.2, while claiming a dramatic evolution in the redshift bin 
1.2-1.5 with blue star-forming galaxies belonging preferen- 
tially to the highest density environments. Since they define 
environmental properties based on a volume limited sample, 
i.e. an absolute magnitude cut at all redshifts, which is bi- 
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Fig. 10. — The U-B rest-frame color as a function of mass, redshift, morphology and environment. The color coding is according to the moiphological 
classes (red-early, green-intermediate, blue-late). The bottom row shows the total values, i.e. the first three rows, from the top, collapsed. Contours define the 
morphological classes and are drawn as explained in the text. The two dotted fines show the estimated mass completeness values of our sample at the upper and 
lower edges of each redshift bin. 



ased against old stellar populations, it would be interesting to 
check whether this effect persists if defining environmental 
densities with a mass complete sample. We postpone a more 
detailed investigation of this issue to a forthcoming paper. As 
of today, only a photometric reshift survey like COSMOS, be- 
cause of the high statistic available, may in fact allow such an 
investigation. 

Cassata et al. (2007) analysed a sample of 2041 galaxies 
with photometric redshifts within 0.61 < z < 0.85 and cover- 
ing an area of 270 arcmin^, that is about 9 times smaller than 
the one we use here. The accuracy of the their photometric 
redshifts is comparable to the one in this work (see Mobasher 
et al. 2007). They perform a morphological classification up 
to Iab < 24, thus the same flux limit we use in our analysis, 
using a set of nonparametric measurements of the galaxy light 
distribution (Conselice 2003, Abraham et al. 2003, Lotz et al. 
2004). Also, following Dressier (1980) they estimate galaxy 
projected densities. They do find a well defined morphology- 



density relation in place and also a clear separation between 
early and late types, in all environments, in the stellar mass 
vs. V-z observed color plane, with bulge-dominated galaxies 
sitting mainly in the red sequence and disk-dominated ones 
populating the blue cloud. Furthermore, they fit and discuss 
extensively the different slopes of the blue cloud (late-type 
galaxies) and the red sequence (early-type galaxies) distribu- 
tions. We see similar characteristics in our data but, in our 
opinion, this is an effect mostly due to the sample mass com- 
pleteness. As it is clear by looking at the mid column of Fig- 
ure 10, the pronounced slope of the late-type galaxies distri- 
bution, compared to the bulge-dominated galaxy population, 
is mostly drawn by the mass completeness: red objects with 
such low mass, if existing, would have an 1 band magnitude 
fainter than 24. 

Nonetheless, their main results are in excellent agreement 
with the one presented in this work at z«0.7, in spite of the 
different ways used to estimate galaxy environmental proper- 
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ties. 

7.4. The specific star formation rate 

One way to explore the contribution of star formation to the 
stellar mass growth in galaxies of different mass, is to study 
the redshift evolution of the specific SFR (SSFR), which is 
defined as the SFR per unit stellar mass. A number of works 
in the last years (e.g. Brinchmann & ElUs 2000, Feulner et al. 
2005ab, Bauer et al. 2005, Iglesias-Paramo et al. 2007, Zheng 
et al. 2007, Cowie & Barger 2008, Damen et al. 2009, Dunne 
et al. 2009, Santini et al. 2009, Pannella et al. 2009) have fo- 
cused on the SSFR redshift evolution, consistently finding an 
increase of the mean SSFRs with redshift. 

Following Madau et al. (1998), the SFR of individual ob- 
jects can be estimated from the rest-frame UV luminosity as 
SFR280() = 1.27 X 10"^*^ X L28()0 in units of M0yr"', where the 
constant factor is computed for a Salpeter IMF. The SFR of 
each galaxy is corrected by the dust attenuation obtained for 
the SED fitting using the extinction curve of Calzetti et al. 
(1994). 

In Figure 1 1 we present the SSFR as a function of stellar 
mass and morphology, in three redshift ranges up to redshift 
1.15, and in different environments. The Figure shows con- 
tour isodensity levels drawn as in Figure 10. 

Before discussing in some detail Figure 11, we would like 
to first discuss some evident similarities with Figure 10. The 
rest-frame U-B color is bracketing the B aimer break of a 
galaxy spectrum, hence is mostly sensitive to the age of the 
galaxy stellar population (e.g. Gallazzi et al. 2005). To con- 
vert an age in a star formation rate, one needs, to first order 
and neglecting the dust attenuation correction, to normalize it 
to the underlying continuum emission or to the galaxy stel- 
lar mass i.e.: LogSFR ex -{U -B)+LogMt,. Taking this into 
account one can, again roughly speaking, link the two plots 
with the following conversion: LogSSFR cx -(U-B). With 
this simple conversion in mind, it is easier to understand why 
there are so many similarities between the two Figures. 

In both Figures the red and blue clouds are very well 
separated, at all redshifts and in all environments. This is 
more pronounced in high density environments, while go- 
ing to lower density regions there is clear trend of the low 
mass early-types to head toward high SSFRs, or bluer U-B 
color, values. These galaxies are young blue bulge-dominated 
galaxies with high star formation activity and stellar masses of 
order a few times lO** solar masses. 

The intermediate-type galaxies fall in between, again con- 
firming their nature of transition objects. It is worth noting 
though that while they overlap with the low mass star forming 
late-type population at aU redshifts and in aU envirormients, 
they tend to leave the high mass end at progressively lower 
redshifts. 

Coming more specifically to Figure 11, the location in the 
SSFR vs stellar mass plot of the maximum galaxy density of 
the early-type population, is increasing its star formation by 
not more than a factor three over the explored redshift range, 
while the blue cloud is moving by approximately a factor ten 
in the same redshift range. The two clouds are pretty much 
segregated in mass: while late-types are always dominating 
at low masses, the early and intermediate type objects com- 
pletely dominate the high mass tail. This mass segregation 
gets definitely stronger in denser environments. 

The upper envelope of the SSFR is running essentially par- 
allel to fines of constant SFR, and it is shifting to higher SFRs 
with increasing redshift, as it was already noted in earlier 



work. Note that this upper envelope is partly affected by a 
selection effect: heavily dust obscured starbursts cannot be 
detected in our sample. 

Furthermore, it is evident that the most massive galax- 
ies have the lowest amount of star formation per unit stellar 
mass up to the highest redshift probed, and hence they can- 
not have formed the bulk of their stars in this redshift range. 
This pushes the bulk of their star formation at earUer cosmic 
epochs, in agreement with the downsizing scenario (Cowie 
et al. 1996). 

Caputi et al. (2008) analysed a sample of 609 sources 
in the COSMOS field, with both optical spectra from the 
ZCOSMOS-DR2 and a MIPS 24^ photometric measurement 
available. In their study they compared different kinds of star 
formation rate indicators, i.e. UV light, line emissions and IR 
luminosity, thus allowing an investigation of the dust extinc- 
tion affecting these diagnostics. As a main result, at least for 
what concerns our work here, they find that the Calzetti et al. 
(1994) attenuation law works well for de-reddening both the 
UV light and the H^, Hq, line fluxes. Although the median Ay 
of our sample, as derived from the SED fitting, turns out to 
be somewhat lower than the one they find for the MIPS 24/x 
selected sample, we note that this might be partially due to a 
selection effect of the Caputi et al. (2008) sample. We post- 
pone a more careful analysis of this issue to a future work, 
which will make use of the recently available Spitzer Space 
Telescope mid-infrared observations (Sanders et al. 2007). 

In Figure 11 it is also worth noting the evolution of the 
SSFR as a function of the environment. At a first glance, 
there is not much difference between the different environ- 
ments. Nonetheless, at all redshifts there is a tendency of the 
early-type sample to present a plume toward the high SSFR 
zone, which is very pronounced at lower environmental den- 
sities. 

8. DISCUSSION AND CONCLUSIONS 

In this work we have studied the evolution of the stellar 
mass content of disk-dominated and bulge-dominated galax- 
ies with respect to the total mass budget up to z ~ 1.2, and its 
dependency from the local environment. 

As expected, we find that early-type objects always dom- 
inate the high-mass end, while at progressively lower stellar 
masses late-type objects increase their contribution to the to- 
tal mass density. 

At each redshift, a transition mass can be identified as the 
stellar mass where the transition from a bulge-dominated to a 
disk-dominated stellar mass budget takes place. Based on the 
good statistics available in this work, we are able to confirm 
that the morphological mix at the high-mass end evolves with 
cosmic time, with the transition mass increasing with redshift. 
At z ~ 0.7 we find a transition mass approximately consis- 
tent with the local value M w 5 x lO'^'M©, while at z 1 the 
late-types and early-types contributions become comparable 
atMw 1 X 10"Mq. 

The stellar mass function of disk-dominated galaxies is con- 
sistent with being constant over the redshift range [0.3-1.2], 
with number densities declining by more than an order of 
magnitude in the mass range Log M, = [10.6 - 12]. On the 
other hand, the stellar mass function of bulge-dominated sys- 
tems shows a pure normalization decHne with redshift. Such 
different behavior of the stellar mass functions of late-type 
and early-type galaxies naturally sets the redshift evolution of 
the transition mass. 

The slope of total number density evolution with redshift 
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Fig. 11. — The SSFR as a function of stellar mass and morphology for three different redshift bins and in different environments. Contours define the 
morphological classes and are drawn as explained in the text. Tilted dotted dashed represent constant star formation rates of 1, 3 and 10 MQyr^' and they are 
supposed to guide the eye. The bottom raw shows the total value. The two vertical dotted lines show the mass completeness values at the upper and lower edges 
of each redshift bin. 



is consistent, within the errors, with being the same in all 
mass ranges, at least up to the common redshift completeness. 
This means that the galaxy stellar mass function, in the red- 
shift range and stellar mass range explored, does not evolve in 
shape but only in normalization with redshift (see also March- 
esini et al. (2008) for similar conclusions at higher redshifts). 
According to our data we do not see any differential evolution 
between the high mass and the low mass end of the galaxy 
stellar mass function. 

We have identified in the photometric redshift distribution a 
large scale structure underdensity at z« 0.5. It gets more and 
more significant with increasing stellar mass. The morpho- 
logical evolution in this underdense redshift bin is delayed at 
all stellar masses: the early-type fraction, both with respect to 
the total number density and to the late-type number density, 
is much lower than those in the contiguous redshift bins, in 
agreement with the expectations of the morphology-density 
relation. We notice that the morphological mix in this under- 



dense redshift bin is, at all masses probed, very close to the 
morphological mix at zwl.l, which has a comparable total 
comoving number density of objects. 

We further explored the redshift evolution of the galaxy 
stellar mass function split by morphology and by environmen- 
tal density. In overdense and medium-dense environments, 
the mass function shape is in good agreement with local de- 
terminations. However, we find that the slope of the mass 
function is anti-correlated with the environmental density: the 
less dense is the environment the more steep is the mass func- 
tion. This trend has also been found locally by Baldry et al. 
(2006). 

In overdense regions, at all redshifts probed, a very massive 
early-type population seems to exist, which overcomes the lo- 
cal ^eW mass function determination: we might be witnessing 
the appearance of giant elliptical galaxies in the highest den- 
sity regions. 

The morphologically-split mass content evolves with both 
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redshift and local density, with a striking feature: at differ- 
ent redshifts, the morphological mix and the transition mass 
appear to mainly depend on the local number density. 

We have presented both the fraction of early-type galaxies 
with respect to the total galaxy population and the ratio be- 
tween early to late-type objects at different redshifts and for 
different environments, finding that a morphology-density re- 
lation is already well in place at redshift 1 . 

We have found that early and late-type galaxies in the stel- 
lar mass vs. SSFR plane (or stellar mass vs. rest-frame U-B 
color) are well separated at all redshifts and in all environ- 
ments. The intermediate type galaxies fall in between, con- 
firming their nature of transition objects. At increasing red- 
shift, the peak of the early-type galaxy distribution in such a 
plane is moving toward higher SSFR by a factor < 3. At the 
same time, the peak of the late-type galaxy distribution shifts 
by approximately a factor ten in the same redshift range. 

The early- and late-type galaxy populations exhibit a sig- 
nificant segregation in mass: late-type galaxies dominate at 
low masses while early-type and intermediate objects domi- 
nate the high mass tail. Therefore, early-type and interme- 
diate galaxies drive the evolution of star formation present 
in massive (> 7 x IO^^Mq) galaxies, while lower mass disk- 
dominated galaxies are mostly responsible for the global de- 
cline of the cosmic star formation from redshift one to the 
local Universe. 

While in general this picture seems to be quite similar in all 
environments, in low density regions there is a population of 
relatively massive, early-type galaxies, having high SSFR and 
blue colors. 

We also explored, with the highly homogeneous dataset 
available, the age of the massive early-type galaxy stellar 
populations as a function of environment. The massive early- 
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type galaxies have similar characteristic ages, colors, and SS- 
FRs, hence a similar formation redshift independently of the 
environment they belong to. The age distributions in different 
environments show a difference in the mean ages (neglecting 
the age-metallicity degeneracy) which is in a good agreement 
with published results. 

Finally, in agreement with previuos work, we have found 
that the galaxies hosting the oldest stars in the Universe pref- 
erentially belong to the highest density regions. 
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APPENDIX 
THE MORPHOLOGICAL ANALYSIS 

In the following we describe the analysis aimed at quantitatively determine the morphological properties of the galaxy sample 
we use in this work. The morphological analysis is based on the ACS imaging of the COSMOS field in the F814W passband, 
and was performed on all galaxies in our sample down to a limit of F814W=24. 

We explain below the details of the adopted procedure and the tests that were performed to assess the robustness of our results. 

Fitting Sersic profiles with GIM2D 

We use the package G1M2D (Simard et al. 1999) to fit PSF convolved Sersic (1968) profiles to the two-dimensional surface 
brightness distribution of each galaxy: 

S(r) = S,e"''[*^''"'"'] (Al) 

where r,, is the effective radius of the galaxy, Sg is the surface brightness at r,., n is the power-law index (Sersic index), and k is 
coupled to n such that half of the total flux is always within r^. For n > 2, k « 2«- 0.331; at low n, Kin) flattens out towards 
and is obtained by interpolation. The adjustable form of the Sersic law has the advantage of parameterizing, through the variable 
exponent n, surface brightness distributions including the exponential radial fall-off of the light profile in bulgeless disks {n = 1), 
and the classical de Vaucouleur profile typical of elliptical galaxies (n = 4). 

During the fitting process, the model is convolved with the user specified PSF, and then compared to the input image. When 
fitting a single Sersic model, G1M2D seeks the best fitting values for the following eight parameters: total flux, half-light radius, 
elUpticity, position angle, center of the galaxy, background level and value of the Sersic index. 

PSF influence on the morphological parameters 

The accurate determination of the image PSF is important for the proper determination of the morphological parameters. The 
PSFs used to convolve the profiles were obtained for each individual ACS tile by stacking about 10 high S/N isolated stars. 
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Fig . A 1 . — Comparison of morphological parameters obtained by fitting the surface brightness of all sources falling in two distinct ACS tiles. For each object, 
the fit was performed on both images, and the resulting Sersic indexes and effective radii are compared in the left and right panels, respectively. The recovered 
median and 16th, 84th percentiles of the A«/n and Af/r distributions are inserted in the two panels. 

Although we did not expect the outcome of single Sersic fits to depend strongly on small variations of the PSF, we carried 
out some tests with varying PSFs on a few COSMOS tiles to assess the influence of the PSF on the retrieved morphological 
parameters. Such tests were performed by using: i) a single universal PSF obtained by stacking all the COSMOS high S/N 
isolated stars; ii) the PSF we derived for the FDF in Pannella et al. (2006), which was observed in the same filter F814W and at 
almost the same depth as COSMOS. We found that the largest differences obtained in the output parameters were of order 3-5 
%, meaning that single Sersic fits are pretty robust against such PSF variations. 

We also tested the accuracy of the PSF itself by fitting all the point-like sources as selected with the SExtractor catalog 
parameters. Specifically, we used the neural network star classifier (CLASS_STAR), the half-light radius, and the FWHM, to 
identify point-like sources. When GlM2D-fitting such point-like sources, the vast majority - but the saturated stars - has a 
best-fit half-light radius smaller than 0.01 pixel, meaning that after the code convolves them with the provided PSF they are 
indeed recovered as point-like sources. 

Results on the overlapping COSMOS ACS tiles 

Since the COSMOS ACS tiles are slightly overlapping, we took advantage of this feature to further test the reliability of our 
fitting procedure as well as to assess many other possible systematics, like for instance PSF variation in time and over the ACS 
field of view, or photometric and astrometric distorsion. 

All sources appearing in two tiles were fitted separately in each tile, and the fit results were compared. The outcome of this test 
is shown in Figure Al, where we plot one against the other the GIM2D output (Sersic index and half-light radius) coming from 
the two fits for all such objects (w 8 000). An overall accuracy of A«/« ^ 0.3 and Ar/r ~ 0.1 is obtained. Taking into account 
that in the tile borders only three out of four dithers are overlapping, and hence that the image rms is greater than in the central 
part of the image, one should consider these uncertainties as safe upper limits. 

Results on simulated objects 

We further tested the accuracy of our analysis by running G1M2D on simulated images. This test was performed by adding 
3000 fake objects, one at a time, to the real image, in a blank, pure-sky region. The simulated galaxies were generated with 
the G1M2D task gsimul, by adding on the selected pure-sky region objects with structural parameters uniformly distributed 
in specified ranges (magnitude in [22.5 - 24.], ellipticity in [0.1 - 0.9], position angle in [0 - 90], effective radius in [0 - 60 
pixel], Sersic index in [0.2-6]). Each single object is first added into the image, then it is extracted with the same SExtractor 
configuration file used for the real images. It is then analyzed from G1M2D exactly as done for the real objects. The input 
magnitudes are generally better recovered for objects with n^er < 2 than for objects with n^er > 2, with a median AM = 0.01 and 
0.06, respectively (see left panel of Figure A2). High Sersic index objects have very extended wings that, depending on the total 
flux, can fall under the sky surface brightness, thus for these objects a lower Sersic index, total flux and effective radii are usually 
recovered (see right panel of Figure A2). We find that the profile parameters of simulated objects are well recovered down to 
F814W=24 (see Figure A3). The recovered ellipticity and position angle, as compared to the input ones, show no trends with 
magnitude and are generally very well recovered (see left panel of Figure A3 for the ellipticity). 

GIM2D vs. GALFIT 

We tested our results by running also the GALFIT (Peng et al. 2002) code on four of the COSMOS ACS patches. GALFIT 
is a 2D galaxy fitting software package more recent than Cj1M2D, and it was designed to extract many structural components 
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Fig. A2. — Left: In the two panels output GIM2D recovered vs. input simulated magnitudes are shown for objects with > 2 (left) and < 2 (right). In 
each panel the median differences with en'ors are shown. Right: GIM2D recovered n,er vs. simulated ones. High Sersic index objects with magnitudes in the 
range 23/24, despite being easily detectable, fall with thek extended wings under the sky surface brightness and suffer from a recovered lower value. 




Fig. A3. — GIM2D test results on simulated objects. The plots show the output vs. input ellipticity (left panel) and /ie (right panel). Points are colored 
according to the simulated objects total magnitudes in red (I<22.5), magenta (22.5<I<23.5) and cyan (23.5<I<24) 

from galaxy images. It is very flexible and able to fit an arbitrary number of galaxies simultaneously on an image, making it an 
ideal tool to fit neighboring objects in a sensitive way and not being biased from external informations such as the SExtractor 
segmentations used by GIM2D. During the fitting process, as for GIM2D, the model is convolved with the specified PSF, and 
then compared to the input image. 

Figs. A4 shows a comparison between the two codes. No systematics are present, the median differences for both the recovered 
Sersic index and half light radius are very close to zero. The scatter is about 20% for the Sersic index and 10% for the half light 
radius. 

One can argue that the COSMOS field is not very crowded and that in more crowded fields one could see some systematic 
differences between the two codes. Still, a similar comparison on deep ACS images of the galaxy cluster A1689 (Halkola et al. 
2006) for a sample of bright (I < 22.5) red-sequence early type galaxies, beside a few obvious cases for which GALFIT performed 
best, showed no significant systematics between the two codes (see also Haussler et al. (2007) for slightly different but still similar 
conclusions on the comparison between the two codes). Of course these considerations apply only to single Sersic fitting, no 
simulations and tests were done concerning bulge/disk decomposition. 
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Fig. A4. — GDVI2D vs GALFIT comparison. In Figure the distribution obtained, on 4 COSMOS tile, of An/n (left panel) and Ar/r (right panel) from the two 
codes. The median is and the width is less than 0.2 for An/n and about 0.1 for Ar/r. 
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Fig. B 1 . — Comparison between the zCOSMOS DR2 spectroscopic redshifts catalog and the photometric redshifts catalog used in this work. The dotted lines 
represent 0.04 and 0.12 rms values, while the soUd line is the bisector, see text for details. 



PHOTOMETRIC REDSHIFTS ACCURACY 

Comparison with zCOSMOS DR2 

In Figure B 1 we show the comparison between the photometric redshifts used in this work and the recently released zCOSMOS 
DR2 spectroscopic redshifts catalog (Lilly et al. 2007). 4506 objects are in common with our catalog. The distribution of 
i^z/{'\-+Zspec) has a median value of 0.01 and an rms of 0.04 with only 2.5% of the distribution having a catastrophic photometric 
redshift, i.e. a Az/(1 +Zspec) greater than 0.2. This nicely confirms the results described in Gabasch et al.(2008), which were 
obtained by using a much smaller spectroscopic redshifts sample. 

Comparison with Ilbert et al., (2009) 

In the left panel of Figure B2 we compare the photometric redshifts used in this work with the ones from the recently released 
Ilbert et al. (2009) COSMOS Legacy catalog for a sample of 40717 objects coimnon to both catalogs. We note that not all 
the objects of our morphological catalog are matched because the Ilbert et al. (2009) sample does not contain the photometric 
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Fig. B2. — Left: Comparison between the photometric redshift catalog used in this work and the one recently published from Ilbert et al. (2009). The dotted 
lines represent 0.05 and 0.15 rms values, while the solid line is the bisector Right: Comparison between the photometric redshift distribution (solid line) used in 
this work (as in Figure 1) and the one (dotted Une) fi'om the Ilbert et al. (2009) catalog for the very same sample. 

redshifts for all the X-ray sources detected in the XMM COSMOS images (Cappelluti et al. 2009). The relative accuracy between 
the two photometric redshift realizations is about 5%, which is in a good agreement with the one expected from the absolute errors 
quoted in the two studies. Finally, only 4% of the common objects have photometric redshifts that differ by more than 20%. 

In the right panel of Figure B2 we show the redshift distribution of our sample obtained using the Ilbert et al. (2009) cata- 
log (dotted line) as compared to the one we are using in the work (solid line, see Figure 1). We note that the two distributions 
agree very well, including the detection of the large scale structure underdensity at z~0.5. 



